Elliptic flow fluctuations in Au+Au collisions at y/s NN = 200 GeV 
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Please note that after these results were reported at Quark Matter 2006 and posted on 
the preprint server it was found that what is reported here as elliptic flow fluctuations, 
should rather be taken as an upper limit on the fluctuations. Further analysis has shown 
that fitting the multiplicity dependence of the q-distribution does not enable one to dis- 
entangle non-flow and fluctuations. The data from the q-vector distrubution does not, 
therefore, exclude the case of zero fluctuations. The remainder of these proceedings we 
CN ■ leave as they were originally reported. 



We report first measurements of elliptic flow (v 2 ) fluctuations for the STAR collaboration 
at middle rapidity in y/s NN = 200 GeV Au+Au collisions. We analyze the multiplicity 
dependence of the flow vector length distribution to disentangle non-flow correlations from 
v 2 fluctuations. We find that the width of the v 2 distribution is approximately 36% of the 
mean v 2 and, within errors, independent of collision centrality. This value coincides with 
eccentricity fluctuations, apparently leaving little room for other sources of fluctuations. 

Introduction: A primary goal of the RHIC program is to collide heavy-ions to create 
a quark-gluon plasma (QGP) and study its properties. Measurements of the azimuthal 
anisotropy of particle production with respect to the reaction plane (i.e. v 2 ) [[1] seem 
to indicate that the matter created in collisions at RHIC behaves as a perfect fluid with 
a viscosity near a conjectured lower bound [E]. This conclusion is based primarily on 
comparisons to hydrodynamic model predictions [HHH]. Uncertainty about the conditions 
at the beginning of the hydrodynamic expansion, however, make the conclusion that the 
matter created at RHIC has a viscosity near the conjectured lower bound ambiguous [ 
3]. Since v 2 reflects the initial spatial eccentricity of the overlap region when two nuclei 
collide, fluctuations of v 2 should reflect fluctuations in the initial eccentricity. As such, 
v 2 fluctuations may provide sensitivity to the initial conditions and help determine the 
viscosity and other properties of the matter created in heavy-ion collisions. 

Analysis method: The n th harmonic reduced flow vector is defined as q n ^ x = -Jj^^i cos(tw/>, 
and q U) y = -y^jSj sin(n0j), where M is the number of tracks and 0j is the azimuth angle 
of a track with respect to the reaction plane [[5]. When enough tracks are used in the 
calculation of q n , the central-limit-theorem (CLT) ensures that its distribution will be a 
two-dimensional Gaussian shifted by \/Mv n in the x direction with widths: 

°l = 2 (l + v 2n ~ 2vl + g n ) , a 2 y = - (1 - v 2n + g n ) ■ (1) 



*For the full list of STAR authors and acknowledgements, see appendix 'Collaborations' of this volume 
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g n represents the broadening of the distribution that arises from non-flow correlations [ 
[5] . The exact direction of the reaction plane is not known, so we calculate the magnitude 
of the flow vector \q n \. If v 2 does not fluctuate, the distribution of \q\ should be given by: 

dN _ i c -i(£±ffi ~ / oj\ k ( q n V l/2fc + l\ I ( q n v n VM \ 
q n dq n ^o x o y A °V \v n VM ) k\ \ 2 / \ ^ / 

where V is the gamma function and I& are modified Bessel's functions. Here we consider 
the n = 2 case. If v 2 fluctuates from event to event, the final \q 2 \ distribution will have 
the form: = / dv 2 f(v 2 )-^^, where f(v 2 ) is the probability for an event to have 

a given v 2 . Zero v 2 fluctuations corresponds to / = 5 (v 2 — (v 2 )) where 5 is the Dirac 
delta function. In our analysis we assume a Gaussian shape for the v 2 distribution: 
/ = , 1 exp (— i^pgg/J j anc i extract the r.m.s. (cr„ 2 ) of the distribution from fits to 

data. Different assumed shapes change the r.m.s. by at most 15% (relative). 
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Figure 1. Simultaneous fits to three dN/dq distributions calculated using the same event 
class but using different fractions of the total multiplicity. The left panel shows a fit to 
simulated data and the right panel shows a fit to preliminary STAR data. 



Non-flow and v 2 fluctuations both broaden the width of the observed q 2 distribution. 
This ambiguity is removed by using different numbers of tracks to calculate q 2 . Since 
the v 2 shift of the \q 2 \ distribution scales with \fM, the broadening of the distribution 
from v 2 fluctuations will scale with \fM. We randomly discard tracks within an event 
class to calculate q 2 with reduced multiplicities. Then we fit the different \q 2 \ distributions 
simultaneously to determine g 2 , a v2 , and (v 2 ). Typical fits are shown in Fig.[T]for simulated 
data (left) and for real data (right). Best fit values agree well with the input parameters 
of the simulation. For our data analysis we use 15 million Au+Au collisions at y/s NN = 
200 GeV measured with the STAR detector [ [6] . 

Eccentricity (standard and participant), impact parameter b, Nb in , N part and multi- 
plicity are calculated from a Monte-Carlo Glauber model [ |7| tuned to reproduce the 
multiplicity distribution of 200 GeV Au+Au collisions. We also use the model to subtract 
impact parameter fluctuations. The a v2 we present represents the width of the v 2 distri- 
bution at fixed b. The correction is carried out according to 8a v2 = a 9v2 5a £vart , where 

OSpart *r 

S<y £ part is extracted from the difference between the width of the e part distribution when 
b is allowed to fluctuate and when it is fixed. 
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Results: The (v 2 ) from our fits are shown versus e part in the top panel of Fig. [2] (left). 
Measuring v 2 fluctuations allows us to remove the major source of systematic uncertainty 
for (v 2 ) [ El E]- {v 2 ) projects to zero for e part = [ [S] giving us confidence that our 
measurements are sensitive to the participant axis (a coordinate system which can fluc- 
tuate to a different angle than the reaction plane [HO]). We also note that v 2 (ft. e part as 
would be expected from hydrodynamic models [[3]. The top panel of Fig. [2] (right) shows 
(v 2 ) versus b compared to a phenomenological model v 2 = 0.034e pari (dN/dy) 1 ^ 3 which is 
closely related to v 2 /e versus ^dN/dy plots [HI]. Our data are consistent with this model. 
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Figure 2. The mean of the v 2 distribution ((v 2 )) (top panels) and the r.m.s. width of 
the the distribution (a v2 ) scaled by the mean (bottom panels). Data are presented versus 
participant eccentricity (left panel) and impact parameter (right panel). Various curves 
are explained in the text. Impact parameter fluctuations have been removed so that the 
data points represent the v 2 fluctuations for a fixed value of b. 



The r.m.s. width over the mean of v 2 (<J v2 / {v 2 )) is shown in the bottom panels of Fig. [2j 
a v2 is corrected for b fluctuations so that the measurements represent the width of the 
distribution for a fixed impact parameter value. a v2 / (v 2 ) is approximately 0.36 ±0.03 and 
within errors independent of centrality. In the top left panel, we also show a comparison 
of the (v 2 ) and a v2 values from this analysis to v 2 from 2- and 4-particle cumulant analyses 
(t>2{2} and t>2{4}) [HIE]. The curves derived from the measured a v2 and (v 2 ) agree with 
the previously measured t>2{2} and t>2{4}. The comparison shows that the differences 
between v 2 {2} and f2{4} can be explained by fluctuations without invoking non-flow 
suggesting that any effects from non-flow on the integrated v 2 are less important. This 
is not necessarily true, however, at higher pt where the ratio of f2{4}/i>2{2} decreases 
with pt- This decrease is likely related to an increase in non-flow correlations from hard 
processes at high px- 

In the bottom right panel of Fig. [2]we compare o v2 j (v 2 ) to nucleon Monte-Carlo Glauber 
calculations of a s / (e) where the eccentricity has been calculated with respect to the reac- 
tion plane (standard) or the participant axis (participant). The e st d calculations can be 
excluded by data, demonstrating that our measurements are sensitive to the initial con- 
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ditions. The relative widths of the t>2 and e part distributions are consistent. This leaves 
little room for other sources of fluctuations and suggests that the conversion of e to t>2 is 
identically efficient for every event. This conclusion is contradicted however by the varia- 

1 /3 

tion of (1)2) /£ with (dN/dy) . This contradiction may indicate that the fluctuations are 
overestimated in the Monte-Carlo Glauber model or that other sources of fluctuation are 
small compared to the eccentricity fluctuations. In principle, different treatments of the 
initial conditions can yield different predictions for the a e j (e). For example, calculations 
of the eccentricity based on Colored Glass Condensate (CGC) yield eccentricity values 
approximately 30% larger than a Glauber model [ 3] . If the geometric fluctuations in the 
CGC calculation are the same as for the Glauber model, then the ratio <j e ,cgc/ '(ecgc) 
will be 30% smaller. This would allow room for other sources of fluctuations, consistent 
with the observation that {vzj/s is not constant. 

Hydrodynamic model calculations with fluctuating initial conditions [d2] are also con- 
sistent with our data when a QGP EOS is used. The use of a resonance gas EOS seems to 
increase the ratio of the rms over the mean above the observed fluctuations. This indicates 
that although we have sought an observable that is sensitive to the initial conditions, we 
may learn about the later evolution as well and may even be sensitivity to the EOS. 

Conclusions: We've disentangled non-flow and fluctuations and extracted a V 2 and (v?) 
versus centrality for Au+Au collisions at 200 GeV. a v ^/ (v 2 ) is 0.36±0.03 and within errors 
independent of centrality. These measurements significantly improve the accuracy of {v^), 
are consistent with 2- and 4-particle cumulant measurements and demonstrate sensitivity 
to the initial conditions. Future work will focus on the system-size dependence, energy 
dependence, and covariance of t> 2 for different sub-events {e.g. partitioned according to 
pseudo-rapidity and p?) ■ 
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